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ABSTRACT 



We present Herschel observations of the isolated, low-mass star-forming Bok globule CB244. It contains two cold sources, a low- 
mass Class protostar and a starless core, which is likely to be prestellar in nature, separated by 90" (~ 18000 AU). The Herschel data 
sample the peak of the Planck spectrum for these sources, and are therefore ideal for dust-temperature and column density modeling. 
With these data and a near-IR extinction map, the MIPS 70 yt/m mosaic, the SCUBA 850 /um map, and the IRAM 1.3 mm map, we 
model the dust-temperature and column density of CB 244 and present the first measured dust-temperature map of an entire star- 
forming molecular cloud. We find that the column-averaged dust-temperature near the protostar is ~ 17.7 K, while for the starless 
core it is ~ 10.6 K, and that the eff'ect of external heating causes the cloud dust-temperature to rise to ~ 17 K where the hydrogen 
column density drops below 10^^ cm~^. The total hydrogen mass of CB 244 (assuming a distance of 200 pc) is 15 + 5 M©. The mass 
of the protostellar core is 1.6 + 0.1 M© and the mass of the starless core is 5 + 2 M©, indicating that ~45% of the mass in the globule 
is participating in the star-formation process. 
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1. Introduction 

Bok globules are nearby, small, relatively isolated molec- 
° ular clouds underg o ing l ow-mass star-formation (e.g ., 
; IClemens & BarvainisI 119881: iLaunhardt & HenningI Il997h . 
. They tend to have only one or two star-forming cores that are 
■ embedded in a larger common cloud. These relatively simple 
' characteristics make these objects ideal for studying the detailed 
processes taking place in low-mass star formation. Specifically, 
the temperature and density structure are fundamental phys- 
ical parameters necessary to understan d core fragment ation 



collap se, and chemical evolutio n (e.g.. IWard-T hompso n et al.l 
; 12007'; 'Stutz et al.' '2008', '2009'; ' Launhardtet al.i i2010). The 
, Herschel (Pilbratt et al. 2010) data cover the wavelength range 
■ which samples the peak of the Planck spectrum for cold sources 
° (6 - 20 K). This wavelength regime is critical for accurate 
modeling of the temperature and densi ty structure in the cold 
environments where stars are born (e.g. JShetty et al.ll20Q9l ). 

The Herschel Guaranteed Time Key Program "Earliest 
Phases of Star-formation" (EPoS; P.I. O. Krause) sample con- 
sists of the Photodetector Array Camera and Spectrometer 
(PACS; ll Poglitsch et"ani201C)|) and the Spectral and Photometric 
Imaging Receiver (SPIRE; lOriflSn et aP 12010^ imaging-mode 
observations of sites of both high- and low-mass star for- 
mation. The low-mass Science Demonstration Phase portion 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



of the sample i s the focus o f this contribution. The source 
CB244 (L12 62: lLvndslll962h is a B ok globule at a distance 
of ~ 200 pc (jHilton & Lahullal [19951) . with an approximate ex- 
tent of ~ 6', or about 0.5 pc. The CB 244 globule contains two 
submm peaks, one associated with a Class protostar located 
at R.A. = 23^25'^46.3^ Decl. = -h74°1739.r', and one associ- 
ated with a starless core located at R.A. = 23^25^^27.1^, Decl. = 
-f74°18'25.3''. The p rotostar drives a molecular outflow (e.g., 
IClemensetaD Il99ll) . while the detection of the CB 244 star- 
less core was first puj p lished as a n additional source inside the 
globule bv lLaunhardi (Il996h andlShirlev etaP (|2QQQ|) and pro- 
duces both an 8 /im (JTobin et al .1120101) and a 24 ^m shadow. 
The YSO and starless core are separated by ~ 90'', and are 
therefore well resolved throughout the PACS and SPIRE bands. 
We use Herschel imaging to construct spatially resolved spectral 
energy distributions (SEDs) of the entire cloud that cover both 
sides of the peak of the SED. These data allow us to reconstruct 
the dust-temperature map and column density distribution of the 
Bok globule and hence the density profiles and mass distribution 
with unprecedented accuracy. They also reveal the role of exter- 
nal heating and shielding by the envelopes in the energy balance 
of such cores in isolated globules. 



2. Observations and data processing 

2.1. Herschel observations 

The source CB 244 was observed with the PACS instrument on 
board the Herschel Space Observatory on 2009, December 30, 
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Fig. 1. 8' X 8' images of CB 244, shown at the wavelengths indicated in the bottom-right of each panel. The images are shown 
on a log scale at stretches intended to highlight the structure of absorption or emission at each wavelength. The locations of the 
protostar and starless core, derived from the 8 yum image, are indicated with x-symbols (note the offset between the x's and the 
starless core emission at 100 - 250 yum wavelengths). The panel contours 'are labeled with their corresponding wavelengths. The 
contour levels are: Ay extinction map levels = {5,10,15,20,25} mags (inner dotted contours indicates the region inside which we 
place lower-limits on the extinction); SPIRE 350 yum = {0.5,1.0,1.5,3.5} mJy/D''; SCUBA 850 yum = {0.3,0.75} mJy/D^'; IRAM 
1.3 mm = {0.1,0.3} mJy/D''. Approximate beam sizes are indicated as yellow circles. 



during the Science Demonstration Program. The globule CB 244 
was observed at 100 yum and 160 yum. We obtained two orthog- 
onal scan maps with scan leg lengths of 9' using a scan speed 
of 20'7s. The scan leg position angles guarantee an almost ho- 
mogeneous coverage of CB 244. We produced highpass-filtered 
intensity maps from these data using the HIPE software package 
dOti 1201 0), version 3.0, build 455. Besides the standard steps 
leading to level- 1 calibrated data, a second-level deglitching 



as well as a correction for offsets in the detector sub-matrices 
were performed. Finally, the data were highpass-filtered, using 
a median window with a width of 27 1 data samples to remove 
the effects of bolometer temperature drifts during the course of 
the data acquisition. Furthermore, we masked emission struc- 
tures, in this case both the YSO and starless core regions, before 
computing the running median. Masking bright sources mini- 
mizes over-subtraction of source emission in the highpass fil- 
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tering step. Finally, the data were projected onto a coordinate 
grid using the photProject routine inside HIPE. As a last step, 
the flux correction factors provided by the PACS ICC team were 
applied. We note that these data have relatively flat background 
values, indicating that the highpass reduction used in conjunc- 
tion with the masking is a relatively robust scheme for avoiding 
artifacts and recovering extended emission in the PACS maps. 
In addition to these processing steps, we checked the pointing in 
the PACS 100 yum map against the MIPS 24 yum mosaic of the 
same region. Using point-sources detected in both images, we 
found a pointing ofl'set of ~ 1.25'' in the 100 jjm map relative to 
the 24 yum mosaic. Because the 160 jim PACS map contains no 
point-sources and was acquired at the same time as the 100 yum 
data, we blindly applied the same pointing correction. The final 
PACS 100 and 160 fim images are shown in Fig.[T] 

Maps at 250, 350, and 500 yum were obtained with SPIRE 
on October 20, 2009. Two 9' scan legs were used to cover the 
source. Two repetitions resulted in 146 s of scanning time with 
the nominal speed of 30'7s. The data were processed within 
HIPE with the standard photometer script up to level 1 . During 
baseline removal, we masked out the high-emission area in the 
center of the field. For these data no cross-scan was obtained; 
therefore, the resulting maps still showed residual stripes along 
the scan direction. We used the Bendo et al. (2010) de-stripping 
scheme to mitigate this eff'ect. In addition, we checked the point- 
ing in the SPIRE data by cross-correlating the images with the 
PACS 160 yum image; we found that the pointing off'sets for all 
three wavelengths are on order of 2'' or smaller, and therefore 
we did not correct for this. The SPIRE 250, 350, and 500 yum 
images are shown in Fig.[T] 



2.2. Other data 

Near-IR data: We observed CB 244 in the near-IR Ks and H 
bands using Omega2000 at the Calar Alto Observatory. The 15'x 
15' field-of-view provided complete coverage of the globule. A 
standard near-IR data reduction scheme was adopted. The ex- 
tracted photometric data where cali brated using t he Two Micron 
All Sky Survey (2MASS; Skrutski e et al.ll2006h . Following the 
NICE method (L adaet al.Ml99 4). the observed (H - Ks) colors 
of stars were then related to visual ext inction by (H - Ks) = 
{(H-Ks)o}-^0.063xAv, assuming the Rieke & Lebof sky| (ll985|) 
extinction law, where ((// - Ks)o) is the mean intrinsic color of 
the stars determined from edges of the field. The Ay values for 
each star were used to create an extinction map. The resulting 
extinction map probes extinction values up to Ay ^ 30 mag, 
with a 3(T error of about 3 mag at Ay = mag. 
Spitzer data: The Spitzer observations presented here are from 
two programs: the MIPS observations are from program 53 (P.I. 
G. Rieke), while the IRAC observations are from program 58 
(P.I. C. Lawrence). The MIPS observations were reduced using 
the Data Analysi s Tool (DAT:lGordon et al.1 ^5) according to 
steps outlined in .Stutz et al.l (l2007l). The IRAC frames were pro- 
cessed using the IRAC Pipeline vl4.0, and mosaics were created 
from the basic calibrated data (BCD) frames using a custom IDL 
program (see Gutermuth et al. 2008). 

(Sub)mm continuum data: The S CUBA 850 ^/m a n d the I RAM 
1.3 mm maps are presented in jLaunhardt et al.l (|201Q|) . The 
850 yum and 1.3 mm maps cover both sources, but do not ex- 
tend very far into the envelope. The maps have native beam 
sizes of 14.9" (850 urn) Sind 1 0.9" (1.3 mm). For further de- 
tails see iLaunhardt et al.l (l2010l) . In Fig. [T] we show the 850 yum 
and 1.3 mm emission as contours. 



3. Dust temperature and column density modeling 

We used the PACS 100 and 160 yum, SPIRE 250, 350, and 
500 yum data in conjunction with the NIR extinction map, MIPS 
70 yum, SCUBA 850 yum, and IRAM 1.3 mm dust emission maps 
to model the line-of-sight-averaged temperature and column 
density of CB 244. The calibrated maps were all homogeneously 
convolved to the largest FWHM beam size-of the data set, the 
SPIRE 500 yum 37'' beam, assuming Gaussian beam profiles, 
and projected onto a common 8''p/x"^ grid. This approach has 
the drawback that all spatial structures smaller than this beam 
size are smoothed, but ensures that the fluxes at all wavelengths 
in a given image pixel are derived from the same area in the 
source. We corrected the filtering of large spatial-scale emission 
in the PACS data using the Schlegel et al. (1998) 100 yum IRAS 
maps and the ISO Serendipity Survey observations at 170 yum: 
we added a zero level of emission of 0.3 and 0.8 mJy/D'' to the 
100 and 160yum maps, respectively. 

For each image pixel, an SED was extracted and fitted 
with a single-temperature modified black-body of the form 
Sy = Q5v(y, Td) (l - e-^^""^), where D is the solid angle of the 
emitting element. By is the Planck function, T^ is the dust- 
temperature, and T(y) is the optical depth at frequency v. In a 
first iteration, we used a simple dust opacity model of the form 
a:^ oc v^ to fit foryS, using ;^^ minimization. We found that j3 =1.9 
- 2 filed the data best over the entire cloud. Note that compact 
structures were smoothed out; therefore this cannot be taken as 
evidence against grain growth in the centers of the two dense 
cores. In a second step, we then fixed the dust model to the 
[Ossenkopf & Henning' ('1994') model for an MRN grain size dis- 
tribution, with thin ice mantles and no coagulation, which has 
a A:i.3mm = 0.51 cm^ g"^ The hydrogen-to-dust mass ratio was 
fixed to the canonical value of 110 (e.g., ?). We then searched 
for the optimum T^ and hydrogen column density A^h values 
(the two free parameters) by calculating flux densities and Ay 
values and then comparing to the emission and extinction ob- 
servations, using x^ minimization. The resulting line-of-sight- 
averaged dust-temperature and column density maps are shown 
in Fig. [21 together with the SED fits at the central positions of 
the two submillimeter peaks. We found the best-fit column- 
averaged dust-temperatures for the protostar and prestellar core 
of ~ 17.7 K and ~ 10.6 K, respectively. Despite the fact that each 
individual image pixel was fitted independently and the maps at 
diff'erent wavelengths have diff'erent outer boundaries, both the 
temperature and the column density maps have very little noise 
and are very smooth, demonstrating the robustness of our fit- 
ting approach. Furthermore, we evaluated the effect of a possi- 
ble filtered-out extended emission component in the chopped 
850 yum and 1.3 mm maps. Adding in off'sets of 60mJy/15'' 
beam at 850yum and 12mJy/ir' beam at 1.3 mm, correspond- 
ing to 20% of the peak surface brightness of the cold prestellar 
core, did not aff'ect the fitted temperature and column density 
(less than a 2% change) for the prestellar core or the protostar. 
Hence, we conclude that possible missing extended emission in 
the (sub)mm maps does not aff'ect our results. 

4. Summary and conclusions 

As shown in Fig. [2l the line-of-sight-averaged dust- 
temperature decreases constantly without any significant jump 
from ~ 17 ± 1 K at the outer boundary of the cloud, where 
A^H drops below 10^^ cm"^ and Ay ^ 0.3 mag, to ~ 13 K in 
the inner cloud at A^h ~ 1 — 2 x 10^^ cm"^ and Ay ^ 3 - 
5 mag. Note that due to the two major drawbacks of our sim- 
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Fig. 2. Dust-temperature (color) and hydrogen column density (white contours) in CB 244. The line-of-sight-averaged mean dust- 
temperature was derived pixel by pixel from modified black-body SED fits to homogeneously beam-smoothed emission maps from 
Herschel PACS and SPIRE, Spitzer MIPS70, and ground-based extinction mapping and submillimeter dust-continuum maps at 0.8 
and 1.3 mm (see text for details). The top color bar indicates the temperature scale. Column density contours are at {0.1 (thick) 0.3, 
0.5, 1 (thick), 2, 3.5, 5, and 7}xl0^^ H/cm^. The peak column density in the prestellar core just reaches 10^^ H/cm^. Two SEDs and 
the respective modified black-body fits are shown for the positions of the Class protostar ("1") and the prestellar core ("2"). The 
yellow dashed ellipses indicate the areas used to calculate source masses (see Sect. 4). The 37'' beam size is indicated by a grey 
circle. White areas represent regions where there is not enough signal-to-noise and/or coverage to get a reliable fit. 



plified approach, beam-smoothing and line-of-sight-averaging, 
the local dust-temperatures at the source centers will be higher 
in the protostar and lower in the prestellar core. Nevertheless, 
the boundaries of the two sources can be derived from the dust- 
temperature map rather than from the smooth column density 
map. We find that T^ = 13 K marks the mean inner dust- 
temperature of the globule, from where the temperature towards 
the protostar rises steeply and that of the prestellar core de- 
creases. We therefore adopt ellipses that follow the 13 K con- 
tour in the dust-temperature map as boundaries for the two 
sources (see Fig. 2). The mean radii of these ellipses are 50'' 
(10,000 AU) for the protostellar core and 76" (15,000 AU) for 
the prestellar core. Integration of the column density map within 
these boundaries yields Mu = 1.6 ± 0.3 M© for the protostel- 
lar core and Mr = 5 + 2 M© for the prestellar core, where the 
uncertainties are derived from an assumed ± 1 K uncertainty in 
the dust-temperature. The total mass of the globule (within the 
A^H = 1 X 10^^ cm"^ contour) is Mr = 15 ± 5 M©. We note that 
Mgas = 1.36Mh. These masses imply that -45% of the mass of 
the globule is participating in the star-formation process. 



Including the new Herschel data, we derive the following pa- 
rameters for the protostar: Lboi ~ 1 .5 L©, Lsubmm/Lboi ~ 4%, Tboi ~ 
62 K, and Menv/Uoi ~ 1-45 M©/L©, (cf. Launhardt et al. 2oS, 
confirming th e Class classification according to IChen et aD 
(1995) and .Andre et al.l (ITOOQ) . Furthermore, the C^^O (2-1) 
FWHM line width for the prestellar core is ~ 0.9 km s"^ (Stutz et 
al., in prep.); using our fitted prestellar core mass and tempera- 
ture, we find that the ratio of gravitational energy to thermal and 
turbulent energy is Egrav / (Etherm + Eturb) - 0-9, a marginally 
sub-critical value which is ex pected for prestellar cores. For 
comparison. iTobin et al.l (1201 Ol) derive a mass of ~3 - 4 M© for 
the prestellar core (integrated over a similar area and scaled to a 
distance of 200 pc) using the 8 jim shadow (see Fig. 1), a good 
agreement given the uncertainties in both mass derivation meth- 
ods. The prominent and extended 3.6 jim coreshine (see Fig. 1), 
originating from dust-grain scattering of the background radi- 
ation field, is an indication of grain growth dSteinacker et al.l 
2010) in the CB244 cloud. These pieces of evidence together 
indicate that the CB 244 globule, and other globules like it, are 
excellent sources in which to study the earliest phases of low- 
mass star-formation. As a next step and in a follow-up paper, we 



A. Stutz et al.: Dust-temperature Structure of CB 244 



will employ 3D-modeling to overcome the effects of line-of- 
sight averaging and beam- smoothing, and to reconstruct the full 
dust-temperature and density structure of CB 244 and the other 
sources in our sample. Measuring reliable temperatures and col- 
umn densities with the Herschel data in a sample of prestellar 
and protostellar cores is a fundamental step towards revealing 
the initial conditions of low-mass star-formation. 
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